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We formulate a semiclassical theory for systems with spin- 
orbit interactions. Using spin coherent states, we start from 
the path integral in an extended phase space, formulate the 
classical dynamics of the coupled orbital and spin degrees of 
freedom, and calculate the ingredients of Gutzwiller's trace 
formula for the density of states. For a two-dimensional quan- 
tum dot with a spin-orbit interaction of Rashba type, we ob- 
tain satisfactory agreement with fully quantum-mechanical 
calculations. The mode-conversion problem, which arose in 
an earlier semiclassical approach, has hereby been overcome. 
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The incorporation of spin degrees of freedom in the 
semiclassical treatment of quantum systems represents a 
considerable challenge to the theorist. The periodic or- 
bit theory (POT), initiated by Gutzwiller through his 
semiclassical trace formula [[jj, has been very successful 
in promoting the research on quantum chaos [Q] and in 
explaining prominent quantum shell effects occurring in 
finite fermion systems in many domains of physics (see, 
e.g., J3|). The growing interest in "spintronics" , i.e., the 
spin-polarized transport and spin dynamics in various 
electronic materials [Q, with the particular scope of de- 
veloping a spin transistor || , make it desirable to formu- 
late also the semiclassical theory including spin. 

Littlejohn and Flynn have extended the POT for sys- 
tems with multi-component wavefunctions and used it for 
the WKB quantization of spherical systems with spin- 
orbit interactions [^) . They did not give an explicit trace 
formula, however, the main obstacle being the problem 
of "mode conversion" which arises when the spin-orbit 
interaction in phase space locally becomes zero. Frisk 
and Guhr have applied the method of Q to deformed 
cavities and promoted the hypothesis that spin-flips oc- 
cur at the mode-conversion points through diabatic tran- 
sitions between the two adiabatic spin-polarized energy 
surfaces. Bolte and Keppeler Q have derived a relativis- 
tic trace formula from the Dirac equation and studied 
nonrelativistic limits, thereby justifying some ad hoc as- 
sumptions made in [Q . Their approach works well in the 
weak-coupling limit where the spin does not affect the or- 
bital motion |J. It may not be extensible |lQ[ , however, 
to strong spin-orbit interactions such as that observed 
in p-InAs or in InGaAs-InAlAs heterostructures JO] , or 
also that amongst the nucleons in atomic nuclei |12| . For 
such systems, the approach presented below provides a 
semiclassical theory with a larger range of validity. 



In this letter, we report on a semiclassical approach 
that allows for the explicit coupling between spin and or- 
bital dynamics. We use the spin coherent state method 
and the path integral in an extended phase space to de- 
rive classical equations of motion which include the spin 
degrees of freedom. We then obtain the amplitudes in 
the semiclassical trace formula for the density of states 
without encountering the mode-conversion problem. 

The path integral for a system with spin in the SU(2) 
spin coherent state representation originally appeared in 
a paper by Klauder [fl3[ as an integral on S 2 . Kuratsuji et 
al. |0 have represented it as an integral over paths in the 
extended complex plane C 1 . From the group theoretical 
point of view, the SU(2) coherent state can be associated 
with the unitary irreducible representation of the SU(2) 
group . To present its construction, we follow a recent 
paper by Kochetov |l6| and his notations. The coherent 
state \z, S) for spin S is defined by 

\z;S) = (l + \z\ 2 )- s exp(zs + )\S,-S), (1) 
S-\S,-S)=0, (2) 

where z S C 1 is a complex number, and s± = si ± i$2, 
and S3 are the generators of the spin SU(2) algebra: 



[h,s±] = ±s± 



[s+,s. 



-] = 2s 3 



(3) 



S € N/2 denotes the representation index. 

The irreducibility as well as the existence of the group 
invariant de Haar measure dfis ensures that the resolu- 
tion of unity holds in the spin coherent state basis: 



z;S)(z;S\dfXs(z) 



'2S+1 



dns{z) 



25+1 



(4) 



tt(1 + |z| 2 ) 2 

which turns out to be the most important property of 
spin coherent states that allows for the path integral con- 
struction, whereby the measure d/xg takes account of the 
curvature of the sphere S 2 . In what follows, we denote 
\z; S) simply by \z) and use z = u — iv with u,v £ R 1 . 

Let us now consider a quantum Hamiltonian with spin- 
orbit interaction 



H = H Q (q. p) + nhs ■ C(q, p) , 



(5) 



where s = (s±, s 2 , S3). Hereby C = {C\,C2,Cz) is a vec- 
tor function of the coordinate and momentum operators 
q, p, and the parameter n regulates the strength of the 
interaction. We make use of Smirnov's idea fl7fl to write 
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the expression for the respective quantum propagator in 
terms of a path integral in both the orbital variables q, p 
and the spin coherent state variables v, u. Imposing pe- 
riodic boundary conditions on the propagation and thus 
integrating over closed paths, we arrive at the expression 
for the partition function (or trace of the propagator): 



Z(T) = J V[q]V[p]V^s[z] exp{z^[q,p,z;T]/^}, 



(6) 



where 1Z is Hamilton's principal action function, calcu- 
lated along closed paths over a time interval T 



R[q,p, z;T\ 



i(pq-qp) 



2Sh 



(uv — vii) 

(i+M 2 ) 



H(q,p,v,u) 



dt, (7) 



and the path integration in (||) is taken over the 2(d + 1)- 
dimensional extended phase space: 

V[q\V\p]V^[z] = Jim n f[ dqj{tk ^ jitk) dMz k ) . 



(l + \z\ 2 ) 2 dH (1- 

— K 



|2\2 



4hS du 

(l + \z\ 2 ) 2 dH 
4hS dv 



dn(v, u) 
du 



C(q,p), 



(l + |z| 2 ) 2 dn(v,u) 
4 dv 



C(q,p), 



. m dHo , fiC , > dc(q, P ) 

Qi = — = — h KhS n(v, u) . 

dpi dpi dpi 



Pi 



&H 

dq t 



dH 



KhS n(v, u) ■ 



9C(q,p) 



(9) 



dqi dqi 

whose solutions define the classical orbits in the extended 
phase space. Note the appearance of h multiplying S ev- 
erywhere in (^|) , which reflects the non-classical nature of 
the spin. The two equations for the spin variables v, u are 
equivalent to h = — n n x C, with the constraint n 2 = 1. 
For spin-boson coupling in the Jaynes-Cummings model, 



equations analogous to ([5|) have been derived in [Q. 

Sugita |l9| has recently given a re-derivation of Gutz- 
willer's trace formula, starting from Z(T) whose Fourier- 
Laplace transform yields the density of states: 



1 r°° 

g (E) = ^ / e lET ' h Z{T) dT . 



(10) 



Following the general arguments of Gutzwiller Q, one 
may evaluate the integrations in (||) using the stationary 
phase approximation, which becomes exact in the classi- 
cal limit 1Z 3> h. The semiclassical approximation of the 
partition function Z(T) then turns into a sum over all 
classical periodic orbits (po) with fixed period T 

z sc (T) = 5> iK -/ fi fv[r,] expi^J^Tyn}, (11) 

po 

where 7Z po are the principal functions (^), evaluated now 
along the periodic orbits in the extended phase space. 
IZpJ [r], T] are the second variations 

rT 



and J — 



-I 



fe=ij=i 

Hereby the time interval T is divided into N time steps 
tji = kT/N, and = z(t&). Note that the periodic 
boundary conditions enable the antisymmetrization of 
the orbital part of the symplectic 1-form in (^) . 

The classical phase-space symbol W(q, p,v,u) of the 
Hamiltonian (^), appearing in the integrand of (Q), is 

W(q, p, v, u) = H (q, p) + kHS n(v, u) ■ C(q, p) , (8) 

where 7io(q, p) and C(q, p) are the Wigner symbols of 
the operators Ho and C, and n = (m, 712,713) = (z\s\z)/S 
is the unit vector of dimensionless classical spin compo- 
nents given in terms of u and v by m — 2u/(l + \z\ 2 ), ^(2) _ 1 
n 2 = 27;/(l + |z| 2 ),andn 3 = -(l-|z| 2 )/(l + |z| 2 ). 2 j-f 

The path integral (|^) receives its largest contributions 
from the neighborhood of the classical paths along which 
the principal function 1Z is stationary according to Hamil- 
ton's variational principle 61Z = 0. The first variation 
hereby yields the following equations of motion 



dt. 



(12) 



where r\ is the following 2(d + l)-dimensional extended 
phase-space vector of small variations 



77 = (A,i/,p,£) 



. . 2Vns . . 2VhS 

oq, ov , op, du ■ 



(i + M 2 )' 



(i + N 2 )i' 



is the 2{d + l)-dimensional unit sym- 



plectic matrix. HS 2 ^ is the second variation of the classi- 
cal Hamiltonian, calculated along the periodic orbits: 
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dqidqj ' *~" lPj ~dqidpj 
^ (-uC 1 -vC 2 + C 3 ) 



d 2 n 



+ > 2 +a 



PiPj 



d 2 n 

dpidpj 



(i + M 2 ) 



<9n ^ dn 

dv du 



dC 



dC 



dp. 



(13) 



After a stationary-phase evaluation of the Fourier integ- 
ral ([lO]), one finally obtains Gutzwiller 's trace formula in 
the standard form p|Jl9) . What is new here is that all 
its ingredients are obtained from the dynamics in the ex- 
tended phase space: the actions S po (E) — lZ po +ET po (at 
fixed energy E) and the periods T po — dS po /dE of the 
periodic orbits, the Maslov indices a po for which Sugita 
has given general formulae p9|| , and the monodromy 
matrix M po defined by rj(T po ) = M po r](0) in terms 
of the solutions of the linearized equations of motion 
i] = J dW^'jdr]. We refer to a forthcoming extended 
paper [Q for the details of our calculations, where we 
also show that in the weak-coupling limit we obtain the 
same results as Bolte and Keppelcr H . 



2 



We shall presently illustrate our method by apply- 
ing it to a simple model Hamiltonian describing a two- 
dimensional electron gas (S = 1/2) in a semiconductor 
heterostructure, laterally confined to a quantum dot, in- 
cluding a spin-orbit interaction of Rashba type pl[ 



H 



(Pl+fy) 



2m* 



+ V(x,y) + kTl (<7 y p x 



cPy) 



(14) 



Here m* is the effective mass of the electron, <ji are the 
Pauli matrices, and the lateral confining potential V(x, y) 
is approximated as a deformed harmonic oscillator 

V(x,y) = m*(uj 2 x x 2 +u ) 2 y y 2 )/2. (15) 

We use oj x = 1.56 lo , oj y = 1.23 lu , and units such that 
h = m* = l)q = 1 and that E and k become dimension- 
less. The spin-orbit coupling strength k in ( [l4| ) depends 
on the band structure p2[] . E.g., for a InGaAs-InAlAs 
quantum dot with ~ 100 confined electrons one would 
obtain a value of k ~ 0.25. We investigate in the fol- 
lowing a situation where k is large enough for the spin 
dynamics to affect the orbital motion. For the examples 
below we have chosen k = 0.67. The equations of motion 
(^) were solved numerically. 

For a large range of parameters with < k 0.75, 
we find the following periodic solutions of Eq. (||): 1) 
two pairs of adiabatic orbits A x and A y , librating along 
the x and y axes with fully polarized spin n y = ±1 and 
n x = ±1, respectively; 2) two pairs of diabatic orbits, 
and D^ 2 , oscillating around A^r with n y ^ 0; and 3) two 
pairs of diabatic orbits, and D^ 2 , oscillating around 
A^ with n x ~ 0. The superscripts (±) of the diabatic 
orbits denote their senses of rotation. The stabilities of 
these orbits will be discussed elsewhere j2(J . For stronger 
couplings k > 0.75 and for energies E < 8, new orbits 
bifurcate from the A and D orbits. In these situations, 
the stability amplitudes will have to be regularized by 
suitable uniform approximations [p3| . 

The diabatic orbits obtained in our approach reflect 
the explicit coupling of the spin and orbital degrees of 
motion. Such orbits have not been discussed in earlier 
semiclassical methods. The adiabatic orbits with "frozen 
spin" correspond to those discussed in (|J^,|l^| where they 
could not be used semiclassically, though, because of the 
mode conversion occurring at their turning points. We 
repeat that we do not encounter the mode-conversion 
problem here, and that the stability amplitudes of all 
orbits can be readily calculated numerically. 

In Fig. [l] we present the (x, y) shapes of the orbits A+, 
D+l, and D+ 2 (left panels), and the time dependence of 
their spin components n x , n y , and n z over one period 
(right panels), all evaluated at E = 60. (The orbits A^ , 
Dy 1; and T) y2 have analogous shapes concentrated along 
the y axis, and the behavior of their n x and n y compo- 
nents is reversed.) We see that along the diabatic orbits 
and D+ 2 , the spin rotates mainly near the (n x ,n z ) 
plane (with n y ~ 0) in a non-uniform way. The diabatic 
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FIG. 1. Periodic orbits in the 2-dimensional harmonic oscil- 
lator with Rashba spin-orbit interaction (see text for param- 
eters). Left panels: orbits in the {x,y) plane. Right panels: 
Spin components n x , n y , and n z versus time. From top to 
bottom: Adiabatic orbit A+ along x axis with polarized spin 
in y direction, and diabatic orbits and D^" 2 oscillating 
around the x axis with spin rotating near the (n x ,n z ) plane. 



spin-flip hypothesis of has thus been replaced here by 
the more sophisticated spin dynamics obtained from the 
coupled equations of motion (^) . 

In Fig. U we show the oscillating part of the density 
of states Sg(E), obtained quantum-mechanically (solid 
line) using exact diagonalization of the Hamiltonian (|T|) , 
and semiclassically (dashed line) using Gutzwiller's trace 
formula. Since the periodic orbit sum generally does 
not converge in systems with mixed classical dynam- 
ics, we have convoluted it with a normalized Gaussian, 
exp{ — (E/j) 2 } /7 \/7r ■ This brings the sum to convergence 
|f24| , whereby only the orbits with shortest periods con- 
tribute, and reflects the prominent gross-shell structure 
of the quantum spectrum. We have used 7 = 0.6 where 
it was sufficient to include the above twelve primitive 
periodic orbits. We observe a rather good agreement be- 
tween the semiclassical and quantum-mechanical results. 
The regular beat-like structure comes about through the 
interference of the six types of classical orbits which all 
have frequencies close to either uj x or uj y . 
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FIG. 2. Oscillating part of the coarse-grained density of states (Gaussian averaging parameter 7 = 0.6) versus energy for the 
same system as in Fig. |l[ Solid line: quantum-mechanical result, dashed line: semiclassical result using the 12 primitive orbits. 

In summary, we have presented a novel approach to include a non-adiabatic coupling of spin and orbital dynamics 
into the semiclassical theory of periodic orbits. Our approach overcomes some of the difficulties with earlier approaches, 
in particular the restriction to purely adiabatic spin motion and the problem of mode conversion. For the simple model 
of a two-dimensional quantum dot with Rashba spin-orbit interaction, we obtain a satisfactory semiclassical description 
of the coarse-grained density of states. To study its experimental implications, we are in the process of calculating 
semiclassically the conductance of this system in response to an external magnetic field pp|] . 

We should point out that there exists a subtle problem connected with the fact that the measure T>[rf\ of the 
path integral over the extended phase-space variations J7 in ( fTl| ) has the proper normalization only in the large-spin 
limit. For small spin, like S — 1/2, this calls for an appropriate renormalization scheme p5| ]. Although the different 
renormalizations needed for pure spatial motion p6[ | and pure spin motion p3| are known, a valid scheme for the 
combined spin-orbit dynamics under investigation here is yet to be found. It might lead to extra phase corrections in 
the semiclassical expressions, as discussed in [ p7[ , and thereby affect the numerical results to some extent. 

We acknowledge helpful and critical discussions with O. Zaitsev and encouraging comments by K. Richter. 
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